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ABSTRACT 

Future extragalactic sub-millimetre and millimetre surveys have the potential to pro- 
vide a sensitive census of the level of obscured star formation in galaxies at all redshifts. 
While in general there is good agreement between the source counts from existing 
SCUBA (850^m) and MAMBO (1.25 mm) surveys of different depths and areas, it 
remains difficult to determine the redshift distribution and bolomctric luminosities of 
the sub-millimctrc and millimetre galaxy population. This is principally due to the 
ambiguity in identifying an individual sub-millimetre source with its optical, IR or ra- 
dio counterpart which, in turn, prevents a confident measurement of the spectroscopic 
redshift. Additionally, the lack of data measuring the rest- frame FIR spectral peak of 
the sub-millimctrc galaxies gives rise to poor constraints on their rest-frame FIR lumi- 
nosities and star formation rates. In this paper we describe Monte-Carlo simulations 
of ground-based, balloon-borne and satellite sub-mm surveys that demonstrate how 
the rest-frame FIR-sub-mm spectral energy distributions (250-850^iin) can be used to 
derive photometric redshifts with an r.m.s accuracy of ±0.4 over the range < z < 6. 
This opportunity to break the redshift deadlock will provide an estimate of the global 
star formation history for luminous optically-obscured galaxies (Lfir > 3 x 10 12 L©) 
with an accuracy of ~ 20 per cent. 

Key words: submillimetre — surveys — cosmology: observations — galaxies: for- 
mation — galaxies: evolution — galaxies: high-redshift 



1 INTRODUCTION 

The star-formation history of the high- z starburst galaxy 
population can be determined from an accurate measure- 
ment of the integral sub-millimetre (sub-mm) and millimetre 
(mm) source-counts, and the luminosity and redshift distri- 
butions of_^iib i min_andjmT>scIc^te 



bution of^h^sc^c^rces^o^h^^o^al^^ 

measur Edby COBE f j Puget et al. 19961: lHauser et al. 1998T 



Fixsen et al. 1998), places an additional strong constraint 
on the acceptable evolutionary models. It is now generally 
believed that the luminosity density, due to star formation 
activity in galaxies, is roughly constant b etween z ~ 1 — 4 . 
A series of ground-based 850/jm surveys ( 5mail et al. 1997 ; 
Hughes et al. 1998|; [Barger et al. 1998|; |Eales et al. 19~ ; 



Eales et al. 200q ; |5cott et al. 2001| ; fox et al.20ul 



Borvs et 



al. 2001 ), undertaken with the SCUBA camera ( Holland et 
al. 1999) on the 15-m JCMT and 1.25 mm surveys ( |Bcrtoldi 
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et al. 2000) using the MAMBO camera (Kreysa et al. 199£ 



on the 30-m IRAM telescope, have contributed significantly 
to this understanding. Prior to the first SCUBA results, the 
optical surveys suggested (contrastingly) that the density of 
star-for mation declined by a factor of ~ 5 over the same 



rays (typically 5sq. arcmin); and low system sensitivity (a 
combination of instrument noise, size of telescope aperture 
and telescope surface accuracy, sky transmission and sky 
noise) which restrict even the wide st and shallowest sub- 



redshift | range flSteidel et al. 199$ |Madau et al. 1996|). This 
discrepancy has highlighted the importance of including a 
correction for dust extinction in the optical-UV estimates 
of star formation rates. Nevertheless, the contribution of 
optically-selected starburst gal axies to the sub-mm back 



mm survey s to areas < 0.1 deg (Scott et al. 2001; Fox 



et al. 2001). In the effort to obtain these wide-area shal- 



ground still rem ains uncertain (Chapman et al. 200C 



Pea- 



cock et al. 200C ). Hereafter, for the sake of brevity, sub-mm 
and m m wavelengths will be referred to as sub-mm unless 



explicit y stated otherwise, and a sub-mm galaxy is consid- 
ered to be an optically-obscured starburst galaxy detected 
in a bl ank-field SCUBA or MAMBO survey. 



low surveys, the current sub-mm and mm observations are 
necessarily only sensitive to the most luminous star-forming 
galaxies (Lpm > 1O 12 L , or SFR > lOOM yr _1 ) assuming 
the population is dominated by galaxies at redshifts > 1. 

The discrepancy shown in Figjl] between the bright-end 
source-c ounts in the 850 urn SCUBA and 1.25 mm MAMBO^ 
surveys (kales ct al. 200d ; |Bertoldi et al. 2000| ; |Scott et al 



2001 



Borys et al. 2001), and a visual inspection of their 



reconstructed maps, suggests that the clusteri ng of sub- 
mm galaxies m ay be influencing the statistics jGaztanaga 



On jc important caveat to this opening statement is fc Hughes 2001] ). Taken at face value, the steepening of the 



the fact that we actually have no accurate information on 
the redshift distribution of the sub-mm galaxy population. 
There are ~ 100 galaxies identified in 850/im surveys at a 
level > 2 mjy. The lack of confident optical, or IR identi- 
fications (e.g. Hughes et al. 1998; Richards 1999), and thus 
precise redshifts for the sub-mm galaxies, was a source for 
the early scepticism of the claims that there existed no evi- 
dence for a decline in density of star formation at z > 2. 

The use of the radio-sub-mm (850/im/1.4GHz) spec- 
tral index has gained popular i ty as a valuable diagnostic 



of reds hift (Hughes et al. 199 
& Yun I200C ; Email et al. 200 



Carilli & Yun 1999 



Carilli 



Bargcr, Cowie fc Richards 



200C |; |E)unne, Clements fc Eales 200C| ). The radio follow-up 
of SCUB A surveys, and mor e recently 1.25mm MAMBO 
surveys flBertoldi et al. 2000| ), with the VLA at 1.4 GHz 
has become the commonly accepted method to measure the 
redshift distribution of sub-mm galaxies. However this tech- 
nique requires extremely sensitive observations if radio coun- 
terparts are to be detected. For example, in order that a 
radio survey detects > 90 per cent of bright sub-mm galax- 
ies (iS850/mi > 8mJy) at z > 1.5 , then a 3cr < 7/iJy s urvey 
must be conducted at 1.4 GHz (Aretxaga et al. 2002). Un- 
fortunately the necessity of such an ultra-deep radio survey 
prohibits this technique from being applied to the wider-area 
(~ 1-10 deg 2 ) sub -mm surveys that will be conducted in the 
next few years (section 2) without an excessive committment 
of radio telescope time (e.g. 400-4000 hours per survey) . An 
alternative, and, we argue, a more reliable method exists: 
namely the use of sub-mm-FIR photometric redshifts, pro- 
vided that one can obtain experimental data with sufficient 
S/N and absolute calibration accuracy. 

Although the possibility of conducting cosmological sur- 
veys at sub-mm and mm wavelengths has been realised in 
the last few years with the development and successful com- 
missioning of sensitive bolometer arrays, the existing 850 /im 
SCUBA and 1.25 mm MAMBO surveys are limited in their 
ability to constrain the evolutionary models of the sub-mm 
galaxy population. The practical re asons for thes e limita- 
tions have been described elsewhere ( Hughes 2000] ) and can 
be summarised as follows: restricted wavelength coverage 
(enforced by the limited number of FIR-mm atmospheric 
windows available to ground-based observatories); low spa- 
tial resolution (resulting in both a high extragalactic con- 
fusion limit and poor positional accuracy); restricted field- 
of-view with the current sub-mm and mm bolometer ar- 



850/im counts for sources > 10 mjy can be interpreted as 
evidence for an under-density of galaxies in the survey-field 
providing those particular data (due to cosmic variance), or 
perhaps a high- luminosity cut-off in the luminosity function 
for galaxies > 10 13 Lq (Fig.jj]). A plausible alternative expla- 
nation is simply the fact that, with so few bright sources de- 
tected (5-10 galaxies in 0.1 deg 2 ), the counts in this regime 
are extremely sensitive to differences between the SEDs of 
the most luminous sub-mm galaxies. Sensitive, wider-area 
(> 1 deg 2 ), rest-frame FIR surveys with accurate redshift 
information (similar to those discussed in this paper) can 
distinguish between these possible alternatives. 



1.1 Ambiguities in the Counterparts of Sub-mm 
Galaxies 

The current 850£tm SCUBA and 1.25 mm MAMBO surveys 
are struggling to identify the sub-mm sources with their op- 
tical, IR and radio counterparts. Positional errors of ~ 2 — 3 
arcsec are associated with the brightest sub-mm sources 
(•S'ssoMm > 8 mjy) which, with some except ions, are detected 
in the widest, sh allowest SCUBA surveys (Scott et al. 2001; 
Eales et al. 200C| ). These errors can b e improved with follow- 
up mm-inte rferometric observations ( grayer et al. 1998 ; |Gear 
et al. 2000; Lutz et al. 2001). However ambiguous optical 



identifications still remain, even with < 2 arcsec resolution 
and sub-arcsec positional errors (e.g. Downes et al. 1999). It 
should be no surprise that sub-mm selected galaxies, includ- 
ing those with mm-interferometric detections, do not always 
have optical counterparts, since high- z galaxies observed in 
the earliest stages of formation may be heavily obscured by 
dust. Indeed, this is the most compelling reason for con- 
ducting the sub-mm surveys in the first instance. Searches 
for counterparts at longer IR wavelengths have proven to 



be more success f ul (gmail et al. 1999 Frayer et al. 200C 



Gear et al. 2000 Lutz et al. 2001), and have shown earlier 



identifications of lower-z bright optical galaxi es to be the 



incor r ect counterpart to the sub-mm galaxies (Smail et al 



199£; Frayer et al. 2000). The natural consequence of these 



ambiguous and potential mis-identifications is an inaccurate 
determination of redshift distribution, luminosities and star- 
formation history of high-z galaxies. 
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Figure 1. Integrated source-counts at sub-mm wavelengths for a model in which the IRAS 60/^m luminosity function is evolved as (1+z) 3 ' 2 
upto 2 = 2, and maintained constant (with no further evolution) for 2 < z < 6. The SED of Arp220 is adopted to represent the whole 
galaxy population. Left panel: The curves represent the predicted source-counts at 850, 500, 350 and 250/im. The measured source-counts 
from SCUBA surveys at 850/im (Scott et al. 2001, and references therein) are shown as filled squares. Additional 850^im source-counts 
from Borys ct al. (2001) arc shown as solid triangles. The flux-densities of the MAMBO 1.25 mm source-counts are scaled upwards by 
a factor of 2.25 to represent the equivalent measurements at 850/im, assuming the MAMBO sources lie at redshifts < z < 6, and are 
shown as crosses (Bcrtoldi et al. 2000). The models also show the effect of a high-luminosity cut-off (at Lfir > 1O 13 L0) in the sub-mm 
population. The shaded areas illustrate the regions of parameter space that need to be searched with future experiments to improve our 
understanding of the evolution of the luminous high-z sub-mm galaxies. Right panel: Integrated source-counts for extragalactic BLAST 
surveys at 250, 350 and 500/im. The dashed-lines show the same evolutionary model described for the left-panel. The histograms illustrate 
the extracted source-counts at each wavelength from two different Monte-Carlo simulations of 1 deg 2 surveys that use a library of 13 
SEDs. These simulations explore evolutionary models with small differences in the strength of evolution and the redshift distribution of 
sources. The details of the Monte-Carlo simulations are summarised in §2.2 




A(yU,m) X( / u.m) 

Figure 2. Left panel: Rest-frame spectral energy distributions (SEDs) of 13 starburst galaxies, ULIRGs and AGN, normalized at 60^tm. 
Lines represent the best-fitting SED models to the data (diamonds) and include contributions from non-thermal synchrotron emission, 
free-free and grey-body thermal emission. The parameters that describe the composite model for each galaxy, previously unpublished 
SED data and references for the remaining SED data are given in Chapin et al. (in preparation); Right panel: Composite SED models 
for the identical 13 template galaxies shown in the left panel, normalized at 60/im and scaled to the flux density of Arp220 at each 
wavelength. The vertically bounded-regions (slanted, crossed and shaded) correspond to the rest-frame wavelength regimes covered by 
observations at 1.4 GHz, 850^tm and in the BLAST filters (500-250/im) respectively for galaxies at z = — 6. This representation of the 
SEDs shows more clearly the wavelength regimes in which the minimum dispersion in galaxy colours can be expected. Hence photometric 
redshifts determined from flux ratios between ~ 100 — 1000/^m are the most independent of assumptions regarding the SED of the 
sub-mm galaxy. 
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1.2 Uncertainties in the SEDs and Redshifts of 
Sub-mm Galaxies 

Despite the divers ity in the properties of sub-mm galaxies 



( [vison ct al, 200C ) , efforts have been made to place them in 
the context of an evolutionary model of galaxy formation. 
The current estimate of the redshift distribution of sub-mm 
galaxies is based on simple comparisons of their IR-radio 



local starburst galaxies and radio-quiet AGN (Hughes et 


al. 199i 


I Smail et al. 1999; Lilly et al. 199S|: 3arilli & Yun 


199S|; |E 


arger, Cowie & Richards 200C 


: 


Carilli & Yun 20001 


Smail et al. 200C ; |Bertoldi et al. 200C 


Dunne, Clements & 


Eales 2 


300|; |Dunlop 2001; Smail et al. 2001; Chapman ct 



al. 20021). The consensus is that the population of sub-mm 
galaxies is distributed between redshifts 2 = 1 — 4, with 
a median z = 2 — 3. The details of this distribution are 
still unknown, and it is common to describe the redshifts of 
individual galaxies as "in the range of with a best guess 
of....". 

This raises an obvious question: What classes of local 
galaxies offer the most useful analogues of the high-z sub- 
mm population? In order to address this we must measure 
the full rest-frame X-ray to radio SEDs of individual high- 
z sub-mm galaxies, measure the dispersion in the shapes 
of their SEDs and their range of luminosities (which may 
depend on redshift), and thus determine the accuracy with 
which the SEDs of high-z sub-mm galaxies can be charac- 
terised by a limited number of local template galaxies. In 
practise, the observed SEDs of high-z sub-mm galaxies are 
restricted to a few detections and limits at 850, 450^m and 
1.4 GHz, with occasional mm interferometric observations 
and a K-band identification. No data exists at rest-frame 
FIR wavelengths with sufficient sensitivity to place addi- 
tional useful constraints on redshifts. 

As mentioned previously, optical and IR follow-up ob- 
servations of those few SCUBA sources for which unam- 
biguous identifications exist have revealed that the sub-mm 
counterparts are often optically- obscured, extremely red ob 



jects (E ROs) with I - K > 6 ( [Smail et al. 199S| ; |Gear et 



al. 200(j; |Lutz et al. 200l| ), although in a few cases the sub 



mm counterparts are blue galaxies, with a weak AGN (e.g. 
SMM02399-0136, Ivison et al. 2001). Also the suggestion 
that SCUBA galaxies are forming stars at a very high-rate 
(> IOOM0 /yr) implies that local dusty ULIRGs may also 
be suitable analogues of high-z sub-mm sources. Unfortu- 
nately the SEDs of this local class of luminous star-forming 
galaxies vary significantl y with increasing FIR luminosity 
(Banders & Mirabel 1996). Furthermore, there are few star- 
burst galaxies at z > 1 for which complete IR-radio SEDs 
exist. 

These initial SCUBA and MAMBO observations have 
demonstrated the necessity for larger-area surveys, and also 
shorter-wavelength (250-500/^m) sub-mm observations. Al- 
though new (and future) ground-based facilities (100-m 
GBT, 50-m LMT, SMA) will offer partial solutions to some 
of these problems over the next 6 years, there will still re- 
main ambiguity in the optical, IR and radio identifications 
of sub-mm galaxies. The consequence is a lack of precision 
in the redshift distribution and rest-frame FIR luminosities, 
and thus an inaccurate constraint on the star formation his- 
tory of the entire sub-mm population. 



Once it begins operation beyond 2008, ALMA will un- 
doubtedly solve many of these difficulties with its power- 
ful combination of receivers operating between 350/xm and 
3 mm, and high (sub-arcsec) spatial resolution. The ex- 
pected launch of the Herschel satellite in 2008 will also 
bring the chance to map large areas of the extragalac- 
tic sky at 250-500/mi with high sensitivity. However, in 
late 2003 onwards, a series of long-duration (~ 15 day) 
sub-orbital balloon flights will be undertaken by BLAST 
- Balloon-borne La rge- Aperture Sub-millimetre Telescope 
( Devlin et al. 2001 ). BLAST has a 2-m primary aperture, 
and is equipped with large-format bolometer cameras op- 
crating at 250, 350 and 500^im which are prototypes of 
the SPIRE focal-plane cameras for the Herschel satellite. 
BLAST will conduct sub-mm surveys over ~ 0.5 — 40deg 2 
(Table 1), and will provide rest-frame FIR-sub-mm data for 
> 5000 galaxies, with a point-source sensitivity 3 times 
greater than is available from even the largest ground-based 
sub-mm telescopes. 

In this paper (Paper I) we describe Monte-Carlo simu- 
lations of sub-mm surveys that demonstrate how the rest- 
frame FIR-sub-mm SEDs (250-500/im) of galaxies can be 
used to measure their photometric redshifts with a conser- 
vative la accuracy of Az ~ ±0.6 over the redshift range 
1 < z < 6. Thus we can finally break the current redshift 
deadlock that limits our understanding of the evolution of 
the sub-mm galaxy population. From such precision in the 
individual galaxy redshifts, we can determine the star for- 
mation history for luminous sub-mm galaxies brighter than 
3 x lO 12 Z/0 with an accuracy of ~ 20 per cent. We consider 
the specific example of future surveys from BLAST. Further- 
more we illustrate how these estimates can be improved with 
sub-mm observations using the more sensitive SPIRE cam- 
era on the Herschel satellite, and additionally with the in- 
clusion of longer- wavelength groun d-based 850Atm data f rom 
SCUBA. In paper II of this series (Aretxaga et al. 2002) wc 
develop this photometric-redshift estimation technique, and 
apply it to the existing mult i- frequency data for the popu- 
lation of ~ 140 sub-mm galaxies detected in SCUBA and 
MAMBO surveys. We include, for the first time, an appro- 
priate treatment of the observational errors at sub-mm to 
radio wavelengths, and provide an accurate measurement of 
the redshift distribution for sub-mm galaxies. In paper III of 
this series (Gaztanaga et al., in preparation) we address the 
clustering of sub-mm galaxies, given their redshift distribu- 
tion derived from these photometric redshifts. 

The following cosmological model is adopted through- 
out the paper: Qa = 0.7, Qm = 0.3, Ho = 67 kms _1 Mpc _1 . 



2 PHOTOMETRIC REDSHIFTS FROM 

WIDE-AREA SUB-MILLIMETRE SURVEYS 

2.1 Designing the Sub-mm Surveys 

To assess the accuracy of the photometric redshifts that can 
be determined from rest-frame FIR-sub-mm data, we gener- 
ate mock catalogues of sub-mm galaxies from Monte-Carlo 
simulations that mimic realistic observing data (§2.2). Thus, 
we must first have prior knowledge of the areas and sensi- 
tivities of the sub-mm surveys that will be conducted with 
BLAST, SPIRE, or any other existing or future sub-mm ex- 
periment, in order to achieve the scientific goals. We have 




Figure 3. Colour-colour (350/im/500/im vs. 250/im/500^tm) vs. redshift distribution for 424 galaxies detected (> 3a) simultaneously by 
BLAST in all 3 niters in a simulated 1 deg 2 survey. Top Panel: Colour distribution of detected galaxies without the inclusion of obser- 
vational (photometric or calibration) errors. Bottom Panel: Colour distribution of detected galaxies with the inclusion of observational 
errors (5 per cent absolute calibration and a random measurement error drawn from a la = 5mjy distribution), which dominate the 
scatter. The black cross marks the position of a 4 X 10 12 Lq galaxy at z = 2.81. The hashed rectangle represents the la error on the 
measured colours. The redshift probability distribution for this galaxy is shown in FigM 
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simulated realistic surveys to satisfy this requirement. These 
simulations, which include an evolving population of sub- 
mm galaxies that are consistent with the observed number- 
counts (Fig.jl]), and a foreground Galactic cirrus component, 
are described in more detail else where ( Hughes fc Ga ztanaga 



200C 



Gaztanaga & Hughes 2001; Chapin et al. 2001). In par- 



ticular these simualtions help understand how the analysis 
of the experimental survey data can be affected by telescope 
resolution, extragalactic and Galactic source-confusion, sur- 
vey sensitivity and noise, and sampling variances due to clus- 
tering and shot-noise. 

Provided the BLAST surveys are conducted towards re- 
gions of low Galactic cirrus (Jiooum < 1.4MJy/sr), then the 
BLAST survey confusion limit will be dominated by extra- 
galactic sources. In the absence of any foreground cirrus, 
the simulation indicates that the extragalactic 3a confusion 
limit of BLAST at 250-500/xm will be ~ 20 - 30 mjy. For 
comparison the tradition al estimate of extragalactic confu- 
sion, 1 source/30 beams (Condon 1974), suggests that con- 
fusion will begin to dominate below 30 mjy. 

Hence, during a single long-duration balloon flight, 
BLAST will conduct a series of simultaneous large-area sur- 
veys (~ 0.5 - 40 deg 2 ) at 250, 350 and 500/xm, down to a 
la sensitivity of 5 mjy, detecting in each survey ~ 1000 lu- 
minous sub-mm galaxies (> 10 12 L o ) between 0.5 < z < 6 
(Table 1). The accuracy with which photometric redshifts 
can be determined from surveys of this magnitude is dis- 
cussed in the following sections. 



2.2 Monte-Carlo Simulations of Sub-mm Galaxy 
Colours 

Measurements of photometric redshifts determined from 
FIR-sub-mm-radio data usually do not include a rigorous 
analys is of the errors and the uncerta i nties in the calcula- 
tions (iHughcs et al 19981; iBlain 199S|; jCarilli fc Yuri 19991 



Carilli fc Yun 2000|; iDunne, Clements fe Eales 20001; |Eafes 



et al. 200C| ; |Barger, Cowie fc Richards 20001 ). The desire to 
assign a degree of confidence (probability) to the photomet- 
ric redshift of any individual sub-mm galaxy motivated the 
development of the Monte-Carlo simulations, taking into ac- 
count the imprecise knowledge of the luminosity function of 
sub-mm galaxies, the dispersion in the luminosity-dependent 
SEDs, absolute calibration and observational phot ometric 



errors (Aretxaga et al. 2001; Axetxaga et al. 2002). From 
these individual redshifts it is then possible to provide a ro- 
bust statistical measure of the redshift distribution for the 
entire sub-mm population. 

In the case of BLAST (and also SPIRE), the power 
of this simple technique to derive redshifts arises from the 
unique ability of these experiments to conduct observa- 
tions at 250-500/im that bracket the ubiquitous rest-frame 
FIR peak (at ~ 60 - 150/xm) in the SEDs of high-redshift 
(1 < 2 < 4) galaxies undergoing a significant amount of star 
formation. 

The adopted recipe that generates mock catalogues of 
galaxies between z = and z ~ 6, for a given wavelength 
and under the chosen cosmological model, is summarised: 

(i) use a single representative SED to calculate the dis- 
tribution of rest-frame FIR luminosities of galaxies, at all 
redshifts, from an evolving 60/^rn local luminosity function, 
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Figure 4. Redshift probability distributions for a selection of 
galaxies detected by BLAST. In all panels the diamonds mark the 
true redshifts of the galaxies. The photometric redshifts (modes 
and 68 per cent confidence intervals) are indicated in the la- 
bels and represented as the shaded areas of the histograms. 
Top panel: 4 X 10 12 Lq galaxy at z = 2.81, detected at 250/jm 
(19 mjy), 350/xm (21 mjy) and 500/mi (26 mjy), slightly below 
the BLAST survey confusion limit. The colours of this galaxy are 
represented as a cross in FigjiJ Middle panel: 6 X 10 12 Lq galaxy 
at z = 2.81, detected at 250/^m (34 mjy), 350/im (56 mjy) and 
500/^m (45mJy) with high S/N, above the BLAST survey confu- 
sion limit. Bottom panel: 4 X lO 12 Z/0 galaxy at z = 4.46, detected 
only at 500/im (22 mjy). 




Figure 5. Monte-Carlo simulations of the colour— colour (350/im/500^tm vs. 250/im/500^tm) vs. redshift distribution for 1380 galaxies 
detected simultaneously by SPIRE at 250, 350 and 500/mi with observational errors (5 per cent absolute calibration and a random 
measurement error of 2.5 mjy) in a 1 deg 2 survey. The black cross and hashed rectangle marks the colours for the identical 4 X 1O 12 L0 
galaxy at z 2.81 shown in FigsJa and W. 



survey area 


la depth 


no. of pixels 


no. of detected galaxies 


no. of > 5cr 


galaxies 


(deg 2 ) 






> 5cr 


> 10cr 


z > 1 


z > 3 


1.0 


5 mjy 


18334 


835 


265 


765 


147 


2.0 


7 mjy 


36668 


1012 


291 


927 


151 


4.0 


10 mjy 


73336 


1100 


294 


988 


117 


9.0 


15 mjy 


165006 


1111 


247 


1023 


129 


36.0 


30 mjy 


660024 


990 


246 


895 


105 



Table 1. Examples of the number of detected galaxies and their redshift distributions in alternative 50 hour 250/xm BLA ST surveys to 
be undertaken during a series of future long duration (~ 10 — 15 day) balloon flights. 



0[L,O] (Saunders et al. 1990), that reproduces the observed 
sub-mm number counts. We adopt a model of pure luminos- 
ity evolution, <f>[L, z] = 4>[L/ L* (z)], where 



L\z) 



(l + zf- 2 L*(0) 
33.6 L*(0) 



for 
for 



z < 2 

2 < z < 6 



(1) 



It has been verified that an alternative choice of SED can 
compensate for small differences in the evolutionary model, 
and hence still reproduce the source-counts (Fig.^) ; 

(ii) randomly assign an SED from a library of 13 tem- 
plate starburst galaxies, ULIRGs and AGN (FigJ^) to 
each of the galaxies in the mock catalogue, and use this 
SED to calculate the intrinsic sub-mm colours. The most 
crucial assumption of the paper occurs at this point in 
the simulation. We assume no correlation between the 



shape of the SED and luminosity (or redshift), and in- 
stead rely on the assumption that our library of 13 tem- 
plate SEDs is representative of those that exist at high- 
z. The library, which includes M82, NGC1614, NGC3227, 
NGC2992, NGC4151, NGC7469, NGC7771, IZwl, Mkn231, 
Mkn273, Arp220, IRAS10214 and the "Cloverleaf" quasar 
(H1413+117), contains galaxies with FIR luminosities in the 
range of log Lpm/ Lq = 9.6 — 12.3. The sample includes 
a mixture of AGN and starburst galaxies, and composite 
galaxies with both an AGN and a luminous extended FIR 
starburst region. The dust temperature (~ 30 — 60 K) and 
emissivity index (f3 = 1.4 ± 0.3), which define the shape of 
their rest-frame FIR-sub-mm spectra are independent of the 
type or strength of nuclear activity. As mentioned in §1, the 
sensitivity of the sub-mm observations restricts the surveys 
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Figure 6. A comparison of the redshift probability distributions for two galaxies observed in a SPIRE survey, and a shallower BLAST 
survey (previously discussed in Fig.0). Symbols and labels are the same as in Fig^. Left panel: 4 X 10 12 Lq galaxy at z = 2.81 detected 
simultaneously with 3 SPIRE filters (shaded histogram), compared with the BLAST redshift distribution (open histogram). Right panel: 
4 X 10 12 Lq galaxy at z = 4.46 detected simultaneously in 2 SPIRE filters; at 350/im (llmjy) and 500^m (20mjy), compared with the 
distribution obtained from a single filter (500/^m) detection with BLAST. 



to a detection threshold of logLpiR/LQ > 12. Thus, an im- 
plicit assumption in this method is that one can simply scale 
the SED of a lower luminosity galaxy to represent that of a 
high-redshift sub-mm galaxy. Furthermore, there is no con- 
clusive evidence that the rest-frame FIR-sub-mm SEDs of 
galaxies vary with luminosity, and redshift beyond the range 
of SEDs considered here - there is simply a lack of sufficient 
high signal-to-noise observational data, over the appropriate 
wavelength regime, to come to any secure conclusions; 

(iii) add random noise (drawn from the distribution de- 
scribing the measurement errors) , and add a constant abso- 
lute flux calibration error to the intrinsic fluxes to allow a 
more 

racy: la photometric errors of 5 and 2.5 mjy for the BLAST 
and Herschel/SPIRE observations respectively, and in both 
cases an absolute calibration error of 5 per cent. The com- 
bined simulated images of extragalactic sources and cirrus, 
described earlier in §2.1, determine the depth of the surveys 
considered here, ensuring that galaxy confusion dominates. 
However this Monte-Carlo treatment of the photometric ac- 
curacy does not consider the contaminating effects of cirrus, 
clustering of galaxies, shot-noise or projection effects. These 
considerations are deferred to a future paper; 

(iv) finally, extract flux-limited catalogues from a 1 deg 2 
confusion-limited surveys at 250, 350 and 500/im with lcr = 
5mJy and 2.5 mjy sensitivities, appropriate for BLAST and 
SPIRE observations. 

Fig.[] shows an example of one of the possible combina- 
tions of sub-mm colours, for a catalogue of ~ 424 galaxies 
detected (> 3a) simultaneously at 250, 350 and 500^tm in 
a 1 deg 2 BLAST survey. As expected the dispersion in sub- 



mm colours is dominated by the observational errors and 
not by the differences between the shapes of the template 
SEDs. This is particularly true when the colour ratios are 
determined over a relatively small wavelength interval, as is 
the case for the 250, 350 and 500/xm observations. However, 
as we extend the wavelength baseline over which colours 
are calculated then the effect of differences in the galaxy 
SEDs on the photometric colours becomes increasingly sig- 
nificant (see Fig.^). An example of this important effect can 
be found with the inclusion of longer wavelength data to 
the BLA ST observations (§2.3.3), or perhaps more ex treme, 
with the com bination of 850/im and 1.4GHz fluxes ( jCarilli 
which are separated in wavelength by a fac- 
tor > 200. Thus, in general, without an adequate treatment 
of observational errors and a representative set of reference 
SEDs, photometric redshifts will be estimated with an over- 
optimistic accuracy. 



2.3 Sub-mm Photometric Redshifts 

Sub-mm photometric redshifts can be determined by cal- 
culating the probability that the the colours (including er- 
rors) of an observed sub-mm galaxy are consistent with the 
colours of every galaxy (with a known redshift) in the mock 
catalogue. The final redshift probability distribution, P[z), 
of any galaxy is then simply the sum of the individual prob- 
abilities from the entire catalogue, or explicitly 

P(z) = a^$( Cl -c ) , (2) 



realistic assessment of the photometric redshift accu- & Yun 2000) 
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Figure 7. Redshift probability distributions demonstrating the 
advantages of follow-up observations that extend the sub-mm 
wavelength coverage. Symbols and labels are the same as in Fig.^ 
Top panel: 6 X 10 12 L Q galaxy atz = 2.81 detected by BLAST in 
all 3 filters (open histogram), previously discussed in Fig.U. The 
redshift estimate improves with follow-up observations at longer 
wavelengths. The shaded histogram illustrates the benefit of de- 
tecting the same galaxy at 850/im (8 mjy) in a shallow sub-mm 
survey (crgso^in = 2.5 mjy), e.g. similar to the UK 8mjy SCUBA 
survey (Scott et al. 2001); Middle panel: This example shows the 



4 X 10 Lq galaxy at 



4.46, initially detected bv BLAST in 



only one filter (500/im), previously discussed in Fig^ (open his- 
togram). Again, the detection at 850/mi (16 mjy) is sufficiently 
discriminating to constrain the distribution to higher redshifts 
(shaded histogram). Bottom Panel: 5 X 10 12 Lq galaxy at = 4.78 
detected at 850^tm (8 mjy) in a shallow SCUBA survey (open his- 
togram). Despite the non-detection of this galaxy by BLAST, in 
all 3 filters (250, 350, 500/im), these shorter wavelength data can 
place an improved constraint on the redshift distribution (shaded 
histogram). 



where a is the normalization constant, such that 
j-^max p^ z ^ z — i wn ere z max = 6. The multi-variate Gaus- 
sian probability distribution, $, for k colours, is given by 



Co 



) = (2n)- k/2 | A" 



1 1 1/2 



exp (-|(ci - c )'A 1 (c i -c )) 



(3) 



where Cj is the colour vector of the ith mock galaxy, such 
that Zi £ [z — dz, z + dz], and Co is the colour vector of 
the test sub-mm galaxy, for which we want to derive its 
redshift distribution. A is the covariance matrix of the colour 
distribution, of elements 



A a/ 3 = ({c a — C a )(C(3 - Cp)) 



(4) 



for two colours a, f3. In the majority of the analysis per- 
formed for this paper, the non-diagonal elements of the co- 
variance matrix are small. When more than two colours are 
considered, for simplicity, the above prescription can be sub- 
stituted with standard x 2 statistics. 

In this analysis, we have adopted redshift bins of width 
0.5 {dz = 0.25). The most probable redshift of the test sub- 
mm galaxy will be given by the mode of the discrete P{z) 
distribution, estimated as the centroid of the Gaussian that 
best fits the 5 bins around the maximum of P(z). The asym- 
metric error bars (z~, z+) correspond to 68 per cent confi- 
dence levels such that P(z)dz = 0.68 and (z+ — z_) is 
minimised. 



2.3.1 BLAST colours 

Fig.^ shows two examples of the redshift distributions for 
sub-mm galaxies detected by BLAST at different levels of 
significance in all 3 filters. A further example of the broader 
distribution derived for a higher- z galaxy detected only at 
500/im is also shown in Fig.W. The accuracy of the redshifts 
determined from BLAST data is described in more detail in 
§2.4, where we also compare the accuracy to those obtained 
from optical-IR colours and the 850pm/ 1.4 GHz spectral in- 
dex. The dispersion in the sub-mm photometric redshifts 
derived from 3 BLAST bands is Az — ±0.6 averaged over 
all redshifts in the range < z < 6. 

2.3.2 SPIRE colours 

An increase in the sensitivity of the sub-mm observations 
(with SPIRE for instance) naturally translates into an in- 
creased accuracy of the redshift distributions. Fig|B| shows 
the increased sample-size (at all redshifts) for the deeper 
SPIRE survey (compared to Fig||). However the greater sen- 
sitivity of the SPIRE observations (la = 2.5 mjy) does not 
significantly improve the accuracy of the redshift distribu- 
tions, compared to those determined from detections in all 
3 BLAST filters. The redshift uncertainties for these sources 
are due to a combination of the absolute calibration errors 
and the dispersion in the template SEDs. However SPIRE 
does provide a far greater redshift accuracy for the sources 
that BLAST detects in only one or two filters with low S/N 
(FigJ^). Thus, if the advantage of the larger aperture of Her- 
schel, the reduced confusion limit and higher sensitivity of 
SPIRE is to be realised, then a clearer understanding of the 
SEDs of high-z galaxies is first required. 
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Figure 8. Left panel: Photometric redshift vs. true redshift relationship for 424 mock galaxies simultaneously detected at 250, 350 and 
500^tm with BLAST in a 1 deg 2 survey. Open symbols show the relationship inferred using only BLAST data to derive the redshifts. Filled 
symbols show the relationships when the redshifts are estimated using colours based on BLAST and complementary 850/xm detections 
from a ground-based survey. The addition of 850/^m measurements significantly increases the accuracy at z > 4 since, at these redshifts, 
BLAST filters sample the rest-frame mid-IR to FIR (~ 35 — 100/^m). The longer wavelength data bracket the rest-frame FIR peak of 
the highest-2 objects, which enhances the diagnostic power for the photometric technique discussed in this paper. Right panel: Example 
of the correction obtained for galaxies at z > 4 when 850/xm observations are included in the photometric redshift analysis. The redshift 
distributions for a 1 X 1O 13 L0 galaxy at z = 4.65 with detections at 250^m (18mjy), 350/^m (29mjy), 500/^m (27mjy) and 850/im 
(20mjy) are shown as open (BLAST detections only) and shaded (BLAST and 850^m detections) histograms. This correction, whilst 
significant, only applies to 4 per cent of the entire sample detected in 3 bands. 



2.3.3 BLAST + 850fj,m colours 

A greater improvement in the accuracy of the sub-mm pho- 
tometric redshifts can be achieved by following-up BLAST 
(or SPIRE) surveys with longer-wavelength survey data (e.g. 
ground-based 850/im SCUBA data), or vice- versa. Figh] 
shows the result of simulations that combine BLAST and 
850^im observations. The SCUBA data include a photomet- 
ric measurement error of 2.5 mj y ; to compare with the shal- 
low UK 8mJy SCUBA survey ( jcott et al. 200 1| ), and an 
absolute calibration error of 9 per cent. The addition of a 
longer-wavelength datum to the colour information has two 
advantages. First, the inclusion of one more colour simply as- 
sists the ability to constrain the redshift. This is particularly 
important for galaxies detected only at 500 and/or 350/im 
(40 per cent of the sources in the BLAST surveys), which 
lie at the highest redshifts. Second, at z > 4 the BLAST fil- 
ters all begin to sample the short-wavelength (rest-frame 
mid-IR) side of the FIR peak. This 850/im extension to 
the BLAST wavelength coverage (250-500/im) ensures that 
there is still a long- wavelength datum that brackets the rest- 
frame FIR peak (at least until z ~ 7). A clear illustration 
of these benefits is shown in Fig.^ where high-z galaxies 
at z > 4, that had previously been confused with lower-z 
galaxies at z — 2 — 3 (using only data at 250-500/im), are 
now symmetrically distributed about the z — z phot regres- 
sion line over the entire range < z < 6, with an improved 
average error of Az ~ ±0.4. This correction applies only to 



a small fraction (4 per cent) of the galaxies detected in all 
3 BLAST bands. This sub-mm method therefore continues 
to provide un-biased estimates of photometric redshifts for 
the most distant galaxies. 

Given this benefit of including longer-wavelength data, 
the deepest BLAST surveys will be selected to overlap the 
largest (0.5-1.0 deg 2 ) surveys carried out by the JCMT at 
850 /im and IRAM at 1.25 mm. The extended wavelength 
coverage (250/im - 1.1mm) will provide accurate redshift 
constraints out z ~ 8. 

The analysis described in this paper simultaneously 
makes use of all of the available multi-frequency data (un- 
like the simple one-colour vs. redshift treatments described 
elsewhere), and takes into account the measurement errors 
and uncertainties in the SEDs of the template galaxies used 
to extrapolate the rest-frame luminosities of the sub-mm 
galaxies from their observed data. 



2.4 Comparison of Optical— IR, Sub-mm and 
Sub-mm— Radio Photometric Redshifts 

The photometric redshifts of galaxies derived from sub-mm 
BLAST data (250, 350 and 500/Ltm), using the method de- 
scribed above, are compared with the true redshifts in our 
mock catalogues. We have demonstrated that an average 
accuracy of Az(submm) ~ 0.6 can be achieved over the 
range < z < 6, although, as Table 2 shows, the scatter in- 
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creases significantly for z > 3.5 for the reasons described in 
§2.3. The photometric redshifts measured from the optical- 
IR SEDs of galaxies in the HDF ( Fcrnandez-Soto, Lanzetta 



& Yahil 1999 



Rowan- Robinson 2001) have an accuracy that 



Figure 9. Top panel: Sub-mm photometric redshifts vs. true red- 
shifts for 189 mock galaxies simultaneously detected at 250, 350, 
500 and 850/im in 1 deg 2 . The overall r.m.s. scatter in the sub-mm 
photometric redshifts, determined from the combined BLAST 
and SCUBA 850/im data, is A zfsubmm) ~ 0.4 in the redshift 
range < z < 6; Middle panel: Sub-mm-radio photometric red- 
shifts vs. true redshifts for galaxies selected from the same cat- 
alogue as in the top panel, using a'po, the 350-1.4 GHz spectral 
index described by Carilli & Yun (2000); Bottom panel: Optical- 
IR photometric redshifts vs. spectroscopic redshifts in the Hubble 
Deep Field - 108 galaxies (Fernandez-Soto, Lanzetta, Yahil 1999 
- diamonds), 156 galaxies (Rowan- Robinson 2001 - crosses). The 
optical-IR photometric redshifts, for galaxies at Ztrue > 2, show 
a dispersion of A z(opt) ~ 0.5 (excluding the outlying sources). 



varies between A z(opt) ~ 0.1 — 0.6 over the same redshift 
range. However we note that the optical galaxies described 
in this HDF analysis are obviously not drawn from the same 
population as the sub-mm galaxies. The accuracy of the op- 
tical photometric redshifts must be much poorer if the tech- 
nique is applied to the optically-obscured SCUBA galaxies. 
The inclusion of longer wavelength data (e.g. 850/im) to the 
BLAST data improves the photometric redshift accuracy 
to Ajz(submm) ~ 0.4, and extends the redshift range over 
which we can measure reliable redshifts to < z < 6 (Fig.^ 
and Table 2). This means that, despite the greater wealth 
of data describing the optical-IR SEDs of galaxies and the 
number of available optical-IR bands, the sub-mm photo- 
metric redshifts of the sub-mm galaxy population can be 
determined with a similar accuracy to the optical-IR pho- 
tometric redshifts of optically-selected galaxies at z > 1.5 
(Fig§. 

We also draw attention to the lowest-redshift bin, 0.5 < 
2 < 1.5, in Table 2, which shows that the addition of 850/im 
data decreases the redshift accuracy, compared to that de- 
rived from the shorter wavelength 250-500/im observations. 
This is due to the fact that the brightness of low-redshift 
galaxies, assuming our evolutionary model for the luminos- 
ity function, will necessarily have low S/N detections in an 
850/im survey with a la sensitivity of 2.5 mjy. Thus the low 
S/N 850/im detections, when combined with the BLAST 
data, allows a greater range of redshifts to be consistent 
with the observed colours. Conversely the addition of shorter 
wavelength SIRTF data with high S /N will improve the ac- 
curacy of our photometric redshift estimates at z < 1.5. 

The redshift accuracies determined from BLAST and 
BLAST+850^j,m data (cols. 3 and 4 in Table §) are compared 
with those obtained from the identical data in our mock 
catalogues using the 350 GHz (850/im) to 1.4 GHz spectral 
index (afH) described by Carilli & Yun (2000). The sim- 
ulations include a photometric error of 7/iJy, and absolute 
calibration error of 3 per cent to compare with the deepest 
VLA surveys. As Figfil clearly illustrates, the photometric 
redshifts obtained with the sub-mm/radio method suffer a 
large scatter (A z(al 5 ®) > 1.5), particularly at z > 2, where 
small differences in a? 5 4 result in large differences in the 
redshift estimates. The inability of the 850-1.4 GHz index 
to determine accurate red shifts at z > 2 has already been 
noted by Carilli & Yun (|2000|). At z < 2.5, photometric 



redshifts using af 4 are still less accurate than those deter- 
mined from the sub-mm (250-850/im) method described in 
this paper. This is due, in part, to the larger scatter of the 
SEDs at radio wavelengths (Fig.^). Naturally, any redshift 
estimated from a single colour index (e.g. af 5 ,?) is also, in 
general, going to a offer a less stringent constraint than the 
use of two or more colours. 



2.5 Measuring the Star Formation History from 
Sub-mm Surveys 

The measurement of redshift probability distributions for a 
large sample of sub-mm galaxies provides the information 
necessary to derive robust estimates of global properties of 
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Figure 10. Left panel: Global star formation rate history of galaxies detected in at least 2 niters from a combined BLAST (3cr = 15 mjy) 
and wide-area SCUBA 850/xm survey (3cr = 8 mjy). Filled diamonds represent the recovered values of the SFR using the photometric 
redshift distributions. We compare these values with the true SFR included in the mock catalogue (shown as a solid line) before 
measurement errors are introduced. The 3a error bars of the reconstructed SFR are dominated by the uncertainty of the SED that 
should be used to derive the 60/im luminosities of the mock galaxies (which have intrinsically different SEDs), and not by the accuracy 
of the redshift determinations. The long-dashed horizontal line shows an estimate of the minimum SFR density for sub-mm galaxies that 
have been detected in the UK 8mjy SCUBA 850/xm survey (Scott et al. 2001), assuming that the SCUBA galaxies lie between 2 < z < 4. 
The short-dashed curves represent the theoretical SFR due to galaxies of FIR luminosities > 10 11— 13 1/0 derived from the evolutionary 
model adopted in the simulations §2.2. Right Panel: Fraction of galaxies, per luminosity interval, detected at all redshifts in the mock 
combined BLAST and wide-area SCUBA survey. The completeness of the survey is ~98 per cent for galaxies > 3 X 10 12 Lq. 
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Table 2. Comparison of the 1<t redshift accuracy, Az, determined 
from a simulated survey with r.m.s sensitivities of 5 mjy at 250, 
350 and 500/im, 2.5 mjy at 850/im and 7^tJy at 1.4GHz using dif- 
ferent photometric colours in redshift bins between 0.5 < z < 5.5: 
Column 2 - optical data (Fernandez-Soto et al. 1999, and Rowan- 
Robinson 2001 in parentheses); column 3 - sub-mm data from 
BLAST at 250, 350 and 500/^m; Column 4 - extended sub-mm 
data from BLAST and additional 850/xm observations; Column 5 
- sub-mm (850/^m, 350 GHz) and radio (1.4 GHz) data. Redshifts 
have been calculated from ct^ 5 ® , denned as the spectral index be- 
tween 350 and 1.4 GHz, using the formalism described by Carilli 
&; Yun (2000). The redshifts in Cols. 3, 4 and 5 have been calcu- 
lated for galaxies drawn from the same mock catalogue. 



the population, such as the SFR history, even if the redshift 
and luminosity of individual galaxies are not known with a 
great degree of precision. The photometric redshift distribu- 
tion of an individual galaxy allows the determination of its 



probable contribution to the total FIR luminosity density in 
a given redshift interval. 

This calculation is performed via 100 Monte-Carlo sam- 
ples drawn from the redshift probability distribution for each 
galaxy in the mock survey. The FIR luminosities of these 
galaxies are then derived from the fluxes in the mock cat- 
alogue and the new synthetic redshifts. The SED from the 
template library that best represents the colours (including 
observational errors) of the whole population is used to cal- 
culate the necessary k-correction. This SED will in general 
be different from that assigned in the original catalogue. 
We then calculate the uncertainty in the k-correction, and 
therefore the error in the FIR luminosity, by randomly se- 
lecting any SED in the template library. Thus we remove any 
bias introduced into the reconstruction of the SFR through 
the choice of a single SED intended to represent the whole 
population of sub-mm galaxies. Finally, the star formation 
density is calculated from the mean of the Monte-Carlo FIR 
luminosities, assuming SFR/A/ey r- 1 = 2.0 x i0- 10 Lftr/L q 
(Hughes, Dunlop, Rawlings 1997). 

A deep BL A ST survey, combined with a wide-area shal- 
low SCUBA survey (Sssopm > 8 mjy), can recover the 
SFR density due to galaxies with Lpir > 3 x 10 12 Lo_at 
0.5 < z < 5.5 with an accuracy of ~ 20 per cent (Fig.|l0[). 
The completeness of the survey is ~ 98 per cent above 
3 x 10 12 L Q (Fig.0). The SFR density derived from the 
mock combined BLAST and SCUBA survey is consistent 
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with th e preliminary re sults from the UK 8 mjy SCUBA 
survey (Scott et al. 2001), and provides a good illustration of 



the advantage of being able to derive photometric redshifts 
for the entire survey sample over a large redshift range. 



3 CONCLUSIONS 

Monte-Carlo simulations of multi-wavelength (250-850/im) 
sub-mm surveys show that the majority of the scatter in 
the derived sub-mm colours will be due, in approximately 
equal parts, to the dispersion in the SEDs of galaxies in a 
library of templates, and the observational and calibration 
errors. Since we can control to some extent the quality of 
the experimental data, it is still the uncertainty in the ob- 
served (or model) SEDs of high-z starburst galaxies that will 
ultimately limit the accuracy of these photometric redshift 
predictions. Theoretical SEDs, based on radiative transfer 



models for high-z starburst ga laxies (Efstathiou, Rowan- 
Robinson & Siebenmorgen 2000), can help in principle. How- 



ever the low S/N, and restricted wavelength coverage of the 
available observational data for high-z starburst galaxies are 
consistent with such a broad range of theoretical models that 
they currently provide little discriminatory power in deter- 
mining the most appropriate SEDs to use as templates. Con- 
sequently, one of the most straight-forward achievements of 
BLAST, but no less important, will be simply to provide an 
accurate empirical model of the temperature-sensitive SEDs 
for 1000's of high-z galaxies at these critical short sub-mm 
wavelengths (250-500/^m). 

The product of the simulations of sub-mm surveys de- 
scribed in this paper is the ability to calculate the proba- 
bility distribution for the redshift of any individual galaxy, 
taking into account observational errors and the uncertainty 
in the appropriate template SED, without the requirement 
to first identify the optical, IR or radio counterpart. These 
simulations demonstrate that the combination of balloon- 
borne (airborne or satellite) short-wavelength sub-mm data 
at 250-500/mi, and longer-wavelength ground-based 850/im 
data, for a statistical sample of galaxies can provide the 
rest-frame FIR luminosity distribution, and hence the star 
formation history of the entire sub-mm population. 

The precision of redshift estimates determined from 
short-wavelength sub-mm observations (250-500^tm) is en- 
couraging, with a la dispersion of A z(submm) = 0.4 — 0.6 
(depending on the availability of longer-wavelength ground- 
based data, A > 850pim) over the range < z < 6. This 
sub-mm precision is comparable to the accuracy of the 
optical- IR photometric estimates of redshift, Az(opt) ~ 0.5 
at z > 2, and significantly better than those provided by 
the 850/mi/1.4 GHz spectral index. An important conse- 
quence of these sub-mm redshift constraints is the ability 
to measure the star-formation of luminous sub-mm galaxies 
(£fir > 3 x 10 12 Z/q) with an unprecedented error of ~ 20 
per cent. 

The advantage of having a constrained redshift proba- 
bility distribution for individual sub-mm galaxies, without 
optical, IR or radio (continuum) counterparts, is that we can 
now determine the likelihood that a redshifted rotational CO 
transition-line falls into the frequency range of any partic- 
ular spectral-line receiver on the next generation of large 
mm-cm telescopes (e.g. 100-m GBT, 50-m LMT). The op- 



portunities to conduct mm-wavelength spectroscopic follow- 
up observations that will provide definitive molecular- line 
redshifts, and dynamical-mass estim ates of sub-mm galax - 
ies, are discussed further in Paper II (Aretxaga et al. 2002). 

The combination of short-wavelength sub-mm data 
from BLAST, SPIRE or other facilities, and large single- 
dish millimetre wavelength telescopes provides a powerful 
combination to break the redshift deadlock that hinders our 
ability to understand the evolution and nature of the sub- 
mm starburst galaxy population. 
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